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Recently we discussed a niultimetric gravity theory containing several copies of standard 
model matter each of which couples to its own metric tensor. This construction contained 
dark matter sectors interacting repulsively with the visible matter sector, and was shown to 
lead to cosmological late time acceleration. In order to test the theory with high-precision 
experiments within the solar system we here construct a simple extension of the parametrized 
^^ \ post-Newtonian (PPN) formalism for multimetric gravitational backgrounds. We show that 



a simplified version of this extended formalism allows the computation of a subset of the 
^^ ■ PPN parameters from the linearized field equations. Applying the simplified formalism we 

\l ' find that the PPN parameters of our theory do not agree with the observed values, but we are 

able to improve the theory so that it becomes consistent with experiments of post- Newtonian 
O , gravity and still features its promising cosmological properties. 

CN ■ I. INTRODUCTION 

> 



n 

^r ' In a recent article [ll] 



—J. _ In a recent article [ll] we presented a multimetric gravity theory containing N > 2 copies 



^^ . . . ^ of standard model matter and a corresponding number of metric tensors g^ ■ ■ ■ g ■ The 
^) \ action was of the type 



N 



S = SG[g\...,g''\ + Y.SMW,'^']. (1) 



7=1 



j^ ■ This structure tells us that the interaction between the different standard model copies with action 

Sm is mediated only by the coupling of the metrics in the gravitational action Sq- Hence particles 
contained in one standard model sector interact with particles from a different sector through an 
extremely small cross section involving the squared Newton's constant. In the astronomical context, 
each type of matter appears dark to observers residing in a different sector. In the Newtonian limit, 
the gravitational interaction was constructed to be attractive within each sector and repulsive of 
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equal strength between different sectors (wliich was shown in [2] not to be possible in the binietric 
case A^ = 2). We derived a simple cosmological model and could show that, due to the repulsion 
between different types of matter, it naturally features an accelerating late-time expansion of the 
universe. 

A wide range of other alternative gravity theories have been proposed in order to model the 
observed late-time cosmology. These include: modified Newtonian dynamics [3] which is not 
a geometric theory; theories which in addition to a metric contain scalar and vector fields to 
describe gravitational backgrounds [J,l5|; ciirvature corrections in metric theories which may appear 



theory 



m 



and area metric gravity 
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or as Ricci scalar corrections in f{R) theories 



3; 



as full Riemann tensor corrections as in 

higher-dimensional models such as [8jj9|; and structural extensions such as non-symmetric gravity 



, \l^. Of course all alternative theories of gravity have to 



be tested against the data available from solar system experiments. For the alternative theories 
mentioned above this issue has been discussed, for example, in |13l4l8l]. 



It is the aim of this article to further test the predictions of our multimetric theory by high- 
precision solar system experiments. For this purpose we will extend the parametrized post- 
Newtonian (PPN) formalism to multimetric gravity and then app ly it to our theory. This elaborate 



formalism was developed mainly by Nordtvedt 
ries, see 



19[ and Will 20| to test single metric gravity theo- 



21[ for a review. It assigns to each gravity theory a set of ten experimentally measurable 
quantities, the so-called PPN parameters /3, 7, ^, ai . . . as, Ci • • • C4- These parameters appear as 
coefficients in a perturbative expansion of the metric tensor and can be computed by a perturba- 
tive solution of the equations of motion. The values of the PPN parameters of a theory are closely 
linked to its physical properties. Most notably, they measure the non-linearity in the Newtonian 
superposition law for gravity, the spatial curvature generated by matter sources, preferred frame 
effects and the failure of conservation of energy, momentum and angular momentum. 

Most PPN parameters have been determined by a wide range of high precision experiments. 
Their values are fixed within very narrow bounds at /3 = 7 = 1 while all other parameters van- 
ish 23] • This means that there is no experimental evidence for peferred frame effects or a fail- 
ure of conservation of energy, momentum or angular momentum. In particular, the bounds are: 
I7 — 1| < 2.3 • 10~^ from Cassini tracking; |/3 — 1| < 3 • 10"'^ from helioseismology; |^| < 10^^ from 
gravimeter data of the Earth tides; |ai| < 10~^ from lunar laser ranging; \a2\ < 4 • 10^'^ from 
the solar alignment with the ecliptic; {a^l < 4 • 10^^'^ from pulsar statistics; IC2I < 4 • 10^^ from 
observations of the binary pulsar PSR 1913+16; IC3I < 10^^ from lunar acceleration; |Ci| < 2- 10~^ 



and I (^2 1 < 6-10 from combinations of the other PPN measurements. Therefore theories for 



which the PPN parameters take significantly different values are experimentally excluded. 

In this article we will show that these high-precision tests can also be used to test multimetric 
gravity theories where the gravitational field is described by several metric tensors each of which 
governs the motion of one type of matter and a corresponding class of observers. In section HIl we 
will construct a simple extension of the parametrized post-Newtonian formalism to multimetric 
gravity. We consider the physically relevant situation of sources formed from a single type of 
matter for which we show that each multimetric gravity theory is characterized by an extended 
set of post-Newtonian parameters, in complete analogy to the standard PPN formalism. We will 
then show in section IIIII that a subset of these PPN parameters can already be obtained from the 
linearized field equations. In section IIVI we finally apply the linearized multimetric formalism to 



our particular gravity model for the accelerating late universe [IJ] and explicitly compute its PPN 
parameters. We find that these do not match the observed values, but we are able to present 
an improved theory from which the correct values are obtained while the cosmological results are 
unchanged. We will conlude with a discussion in section |Vl The appendix contains some important 
details of the computation. 

II. MULTIMETRIC EXTENSION OF THE PPN FORMALISM 

The PPN formalism was originally developed for single metric gravity theories, see the review 



m 



2ll ]. In this section we present a simple extension of this formalism for multimetric gravity 
theories. This extension is contructed to describe the physical situation of the solar system for 
which we will argue that two of the N metrics suffice. These two metric tensors will be expressed 
in terms of the PPN potentials and an extended set of PPN parameters. Comparison of these with 
the standard PPN parameters then enables tests of multimetric gravity by high-precision data. 

The class (fTl) of multimetric gravity theories we consider is restricted by the following four 

riQ 

assumptions which were motivated in full detail in [1|, I4]: 

(i) The field equations are obtained by variation with respect to the metrics 5^^^ . . . ^f^, and so 
are a set of symmetric two-tensor equations of the form Kab = 8TTG]\fTab- 

(ii) The geometry tensor Kab contains at most second derivatives of the metric, which can be 
achieved by a suitable choice of the gravitational action. 

(in) The field equations are symmetric with respect to arbitrary permutations of the sectors 
(g^,^^), which can be understood as a generalized Copernican principle. 



(iv) The vacuum solution is given by a set of flat metrics g^^ = r]ab- (Poincare symmetry for all 
metrics simultaneously implies 5^^ = X rjab for constants A that, invoking the Copernican 
principle for the vacuum, should be equal and can be set to A^ = 1.) 

Cosmological constants are excluded because we are interested in multimetric gravity theories in 
which the accelerating universe is modelled by a repulsive interaction between different standard 
model copies, as we have seen in [l|]. 

Another consequence of the repulsion of matter from different sectors is their separation as the 
universe evolves. Hence we may safely assume that the gravitational field in regions like our solar 
system is dominated by a single type of matter; we formulate: 

(v) Regions exist where the gravitational field is generated by matter sources from a single sector. 



Combining this assumption with the symmetry assumption Imj) guarantees the existence of solu- 
tions in which the metric and the energy-momentum tensor corresponding to the dominant matter 
source are distinct, while all other metric tensors are equal and their energy- momentum tensors 
vanish. We assume that this simplest solution is actually realized, since, as observers in a distinct 
region, we have no further detailed access to the physics of the other sectors: 

(vi) The metric tensors of all other sectors besides the one distinguished by the dominant matter 
source are equal. 

In the following we will indicate all quantities within the distinct sector / = 1 by a superscript '+', 
and all quantities within the other sectors by a superscript '— '. Hence gf^^ = 5^^ and TJ^ = T^^^ 
while g^f^ = g^^ = . . . = g^^ and T^ = T^^ = . . . = T^ = 0. We will now extend the standard single 
metric PPN formalism to this physical situation. 

Basic ingredient of the PPN formalism is an expansion of the geometric background in orders of 
the velocity of the source matter. Using assumption i fz^) this is a weak field approximation around 
the flat vacuum metric 7] in Cartesian coordinates (x'^,x"), 

± _ , , ± _ , ^(l)± , ;^(2)± , (3)± , {4)± . . 

Higher than fourth velocity order 0(4) is not considered. It turns out that not all metric per- 
turbations h^l ~ 0{i) are relevant to describe the motion of test bodies. Moreover certain 
components vanish due to Newtonian energy conservation or time reversal symmetry. We now 
only list the relevant, non-vanishing components of the metric perturbations. These are written 



in terms of the so-called PPN potentials U,Uaf3,Va,Wa,^w,^i ■ ■ -^^^-^t^ and constant PPN 
parameters a^ ,^^,6^,ajj_, fi^,^^,4>i . . .(j)^,pL^ ,v^ as 

hfil^ = 2a^U, (3a) 

h^T = 27^^'^a;3 + 2e^U^p , (3b) 

42^ = ^+ (K. + W«) + a± (y„ - VF„) , (3c) 

/iS^ = -2/3±t/2 - 2t^w + 20f $1 + 20±$2 + 2(^^$3 + 20f $4 + 2^^^ + 2z.±e . (3d) 

The spacetime dependent PPN potentials are Poisson-like integrals over the source matter energy 
density /?, velocity v°^, internal energy pYi and pressure p. Due to the virial theorem, the energy 
density is associated a velocity order /> ~ 0{2) while one assigns pTi^p ~ 0(4). The PPN potentials 
at second velocity order are the standard Newtonian potential U and a tensor potential Uap, 

[/(x°,x) = / dV 7_, '' , f/a/3(x°,x) = / dV — ^-3 ^. (4) 



Similar integrals define the third order vector potentials 

,0 ^/^„ (^0 ^/\ 



F„(.o,x) = /dV^i^^ML^ 



(5a) 

|X — Jy I 

,,, , -^ / j3 ,/5(2;°,x')t'/3(x°,f')(a;a-2;'„)(X;3-4) .^^, 

Wa(x",f) = / d^'x' -z — ^77^ — (5b) 



2l|] for full detail. 



and the fourth order scalar potentials ^Wi *^ii ^2i ^3, ^4, -^i B, see 

The metric ansatz presented above fully reduces to the PPN formalism in single metric theories. 
To see this, one simply drops all superscripts 'it'. Moreover, it is then conventional to choose a 
gauge so that the parameters 9 and u vanish, thus removing the potentials C/q^ and B. Also, 
one may choose a = 1 by absorbing its value into the definition of the gravitational constant. 
In the extended PPN formalism discussed here we have the same amount of freedom. However, 
gauge choices result from diffeomorphism invariance and affect all metrics simultaneously, and 
also a rescaling of the gravitational constant in one sector will affect the gravitational interaction 
between all sectors. Hence only one of the two metrics g^ can be simplified. 

Since we wish to compare the multimetric PPN parameters to experimental data, we turn our 
focus to observers that reside within the distinct sector, i.e., to observers for whom the dominating 
matter source is visible. These observers as well as their visible type of matter are affected by the 
metric g~^ only. Thus g'^ corresponds to the single metric g in the standard PPN formalism. We 
will therefore choose a gauge in which h^ has the standard PPN form, i.e., 6^ = v^ = 0, and fix 
the gravitational constant so that a~^ = 1. The remaining ten PPN parameters contained in h~^ are 



2ll | and the notation we use in the 



then identified with the PPN parameters known from the single metric case by direct comparison. 
The conversion between the standard PPN parameters used in 
metric ansatz ([3]) is given by the relations 

(/.^ = l + 37-2/3 + C2 + e, </>^ = l + C3, </>! = 37 + 3C4 - 2^ , (6) 

/i+ = ^(2^-Ci), a+ = -^(4 + 47 + ai), al = -^(2 + 47 + ai-2a2 + 2Ci-40. 

With this identification we can convert the experimentally measured values of the standard PPN 
parameters to our notation and obtain 

/3+ = 7+ = 0+ = l, ^+ = ^+ = 0, ,^+ = </,+ = 2, 0+ = 3, a+ = -2, ^+ = -^- (7) 

As we already argued in the introduction, these values are fixed by numerous experiments within 
very narrow bounds. Gravity theories with significantly different PPN parameter values are there- 
fore experimentally excluded. 

The values of the extended PPN parameters can now be computed in complete analogy to the 
standard PPN formalism. We will examine this procedure in the following section Hill Conveniently, 
it will turn out that some of the extended PPN parameters can already be obtained from the 
linearized equations of motion. 

III. LINEARIZED PPN FORMALISM 

In the preceding section we have discussed an extension of the PPN framework to multimetric 
gravity theories, which allows tests of these theories by means of available experimental data from 
solar system and astronomical experiments. We have shown that, in addition to the parameters 
obtained from the standard PPN framework, we obtain further parameters that characterize the 
influence of matter sources from one sector on the metric tensors of the other sectors. In this 
section, we will show that a number of the extended PPN parameters can be computed already 
from the linearized equations of motion. This covers the PPN velocity orders up to 0(3) while all 
terms of 0(4) and higher are neglected. We will give a step-by-step recipe for this calculation. 

A. Geometry and matter content 

The starting point for our computation is the most general linearized field equations compatible 
with our assumptions stated at the beginning of section [III Using the convention Gn = 1 for the 



normalization of Newton's constant, these equations read 

K_ab = SttT,, (8) 

with the hnearized geometry tensor 

Kab = E- dPd^ah)p + £ • ^Kb + E • dadbh + M • dpd.hP^'rjab + K • Ohrjab + 0{h^) (9) 

and constant parameter matrices P, Q, R, M, N. Our assumptions pose two restrictions on these 
matrices. First, the hnearized equations of motion ^ are invariant under a permutation of the 
sectors by assumption liii)) . Hence the parameter matrices wiU have the form 

O^-^ = 0~ + (0+ - 0-)6^-^ (10) 

with diagonal entries O^ and off-diagonal entries O^ for O = P,Q,R,M,N. This leaves us with 
a set of ten parameters determined by the underlying (nonlinear) gravity theory. Second, the 
equations of motion are tensor equations by assumption ([|), and so the linearized equations should 



be gauge-invariant. Using the formalism of gauge-invariant perturbation theory detailed in 
finds the invariance conditions 



, one 



(P + 2£) • 1 = (P + 2^) • 1 = (M + £) • 1 = 0, (11) 

where 1 = (1, . . . , 1)* and = (0, . . . , 0)* denote A^-component vectors. Using ()10p . these conditions 
can be written in the form 

P+ + (iV-l)p- + 2Q+ + 2(iV-l)Q- =0, (12a) 

P+ + (iV-l)P"+2P+ + 2(iV-l)p- =0, (12b) 

M+ + {N-1)M~ +N+ + {N -1)N- =0. (12c) 

These equations fix three of the ten parameters in the parameter matrices so that the most general 
linearized curvature tensor consistent with our assumptions is completely determined by a set of 
seven parameters. 

We will now turn our attention from the geometry side to the matter side of the equations of 
motion. Recall that, according to assumption ^, we consider only solutions of the field equations 
in which the gravitational field is generated by matter sources within a single sector, i.e., by a 
single energy- momentum tensor T^ , while all other energy- momentum tensors T^ must vanish. In 
order to solve the linearized field equations dH), we need to expand r+ up to the required order of 
perturbation theory, i.e., to velocity order 0(3). We will use the ansatz 

^o'o=P> T^a = -PVa, r+3 = 0, (13) 



corresponding to a perfect fluid of density p ~ 0{2) and velocity v"' ~ 0{\). 

B. Computation of the extended PPN parameters 

We will now explicitly solve the equations of motion. Omitting all terms in the PPN metric ^ 
corresponding to perturbations of velocity order O (4) we may use the simplified ansatz 

hoQ = -a^^X, (14a) 

h^^ = aixt + atX-, (14b) 

h% = 29^d^dpx - (7^ + ^^)Ax<5a/3 . (14c) 

These expressions are rewritten in terms of the so-called superpotential 

x(x°,x) = - / d?x'p{x^,x')\x - x'\ (15) 

using the relations 

U = -2^X , Ua/3 = dadpx - 2 Ax<^a/3 • (16) 

We have furthermore introduced the notation X^ = Vai: Wa for the vector potentials. The advan- 
tage in using X^ instead of Va and Wa results from the fact that X~ = dadoX is a pure divergence 
and X^ is a divergence-free vector, d'^X^ = 0. These relations follow from the Newtonian continu- 
ity equation dQp + da{pv°') = and the definitions (0), and will be used repeatedly in the following 
computation. 

We begin by performing a (1 + 3)-split of the equations of motion ([8]). Using the energy- 
momentum tensor ansatz (jl3|) we obtain the equations 

K+=8np, ^00 = 0, (17a) 

K^^ = -87TpVa, ^Oa = 0, (17b) 

^i = 0, K-^ = 0. (17c) 

In order to solve these equations, we expand the geometry tensor K^^ given in ([9]) using the PPN 
metric ()14p . In this calculation we once again drop all terms of velocity order 0(4), taking care 
of the fact that time derivatives count as do ~ 0{1). Up to the required order 0(3) the geometry 
tensor then takes the form 

K^o = cfAAx, (18a) 

Ko± = 4AX+ + cfAX^ , (18b) 

iC± = cfAAx5ap + c^Adadfsx , (18c) 



where the coefficients c^ , . . . , c^ are constants which depend hnearly both on the PPN parameters 
and the components of the parameter matrices (jlOp . For a detailed expansion of these coefficients, 
see appendix \M 

We win now determine the coefficients Cj^ , . . . , Cg such that the equations of motion are satisfied 
for arbitrary matter distributions p and v". First, we solve the scalar equations ()17ap . Using the 
relation 

AAx = -2AU = 8irp, (19) 

one can see that these are solved if, and only if, the corresponding coefficients take the values 

4 = 1, c^; = . (20) 

We continue with the tensor equations (|17c|) . Note that AAxda/s is a pure trace term, while 
Adadpx decomposes into a pure trace and a traceless part, 

Adadpx = AA„/3X + \^^X5aP , (21) 

using the traceless second derivative Aq^ = dadjs — 5ai3A/3. In order for the tensor equations to 
be satisfied, both the trace and the traceless part, and thus the coefficients of both terms, must 
vanish, 

cf = cf = 0. (22) 

Taking a closer look at the expansion of the coefficients displayed in appendix [XJ one finds 
that the coefficients 0^,0^, Cg only depend on the PPN parameters a , 7 , 9 . We therefore have 
obtained the six equations ()20p . ()22p for six of the PPN parameters. However, due to the gauge 
invariance conditions l\l'2\] . these are linearly dependent. In order to solve the equations, one needs 
to (partially) fix a gauge by fixing the value of one of the parameters. The standard PPN gauge 
corresponds to the simple choice 9~^ = 0. 

We now turn our attention to the vector equations (|17bp . The equation for A'J^ is easily solved 
using the fact that both the divergence-free and the total derivative part of the curvature tensor, 
and thus both coefficients must vanish, 

C2" = cj = . (23) 

Finally, we consider the equation for Kq^ which yields 

C+AX+ + c+ AX„ = -Snpva = 2AF, = A(X+ + X") (24) 
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using the definition of Va in the second equahty. The equation above is now spHt into pure 
divergence and divergence-free vector terms which decouple and have to be solved independently. 
This results in 

4 = 4 = 1. (25) 

Another close look at the newly obtained equations for Cg and the expansions given in ap- 
pendix |A] reveals that all terms containing a_ drop out due to the gauge invariance conditions ()12p . 
The equations for Cg then turn out to be linearly dependent on the equations we have obtained from 
the scalar and tensor components of the equations of motion, and thus they are solved identically. 
The remaining two equations for Cg can finally be used to solve for the PPN parameters a_^_ . 

To summarize, the linearized field equations in our multimetric PPN framework already are 
strong enough to determine the eight extended PPN parameters a^, 7^, 9^,a_^_. These parameters 
are the solutions of equations (j20p . (j22p . (j23p and ()25p . Given any particular multimetric theory 
consistent with our assumptions, one may use this result as a quick test of solar system consistency, 
simply by comparing the predicted PPN parameters with the experimentally favoured results ([7|). 
Before analyzing the particular theory proposed in [l| in the following section we remark that 
a calculation of the remaining PPN parameters in the extended multimetric formalism requires 
higher order perturbation theory that also covers velocity orders 0(4). In practice this is a very 
lengthy calculation that we do not wish to enter in this article, but it poses no difficulty in principle. 



IV. APPLICATION TO OUR REPULSIVE GRAVITY MODEL 

We will now determine the PPN parameters a^, 7^, 9^,a_^_ of the repulsive gravity model pro- 
posed in [1] by applying the multimetric PPN formalism developed in the previous sections. For 
this purpose we first derive the linearized field equations and determine the parameter matri- 
ces P,Q,RjM.ji!L that appear in the linearized curvature tensor Q. Second, we follow the steps 
detailed in section IIIIBI in order to compute the PPN parameters of our theory explicitly. It will 
turn out that these do not agree with the values obtained from experiments. But this problem can 
be solved, as we will finally show, by simple correction terms that improve the action of our theory. 
The calculated PPN parameters of the improved theory are now consistent with experiment, and 
the theory features the same accelerating late-time cosmology. 
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A. Linearized field equations 

The action of our repulsive multimetric gravity model []| is of the form ([T|) . The gravitational 
part of the action is explicitly given by 

SG[g\ • • • ' 5^] = Y^ y d'^V9^ E (^ + y^'V'Rfj (26) 

with density go = Y[i=i (d^) ) ^^'^ the matter action is a sum of standard model actions 

SM[g', ^'] = I d^xy^LMig', ^'] . (27) 

In the case A^ = 1 when only a single metric is present, the action becomes identical to the Einstein- 
Hilbert action for an appropriate choice of the constant parameters x, y. The field equations are 
obtained by variation and take the form 

kL = 8^T^b (28) 

with the geometry tensor 

I — r 1 ^ ^ 

kL = ^Jgo/g' -^gL E ix + yS''')g''^R''^, +Y.{x + y6''^)R\, 

L J,K=1 J=l 

N 

E 

j=i 



K9m% - aLsy^j) - 9''%agi)j) E(^ + y^'') 2g'p^s''\^,s',^ (29) 

j=i ^ 

+ Ig'^'S'.S'd + \g'''Vis'd + Vis'-'^dpg'''' + S'-'P^^S'-'^apg''' + S'-'\gS''^dpg'''A 



where the connection difference tensors are defined as 

N 



s 



IJi 



jk 



p/i pjj qIJ oljk bJi \ ^ olJi bJ oJfc ('1(W 

jk ^ jk J >-> j ~ *-" jfc ) *-' jk — "77 / _, >-' jk : *-" j ~ '-' jk ■ l"-*"/ 



1=1 



We now derive the linearized field equations using the perturbative ansatz g^^f^ = rjab + /i^ft- 
Note that the connection differences S ik and the Ricci tensor R- • are of first order in the 
metric perturbations h^ , so any terms containing products of two connection differences drop out; 
covariant derivatives acting on connection differences are replaced by ordinary partial derivatives; 
the metric tensors g are replaced by the fiat metric 77 whenever they appears in a product with 
a connection difference or Ricci tensor. These handy rules significantly simplify the computation 
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and one finally obtains the linearized geometry tensor 

N 



kL = E 



J=l 



{2x - {Nx - y)6'^) dPd^X)p + (-^ +\{Nx- y)5'A Uh 



+ [^5'' - - - ^) ^«^^^' + (f^^" - - - ^) d.^h^''^'^^ (31) 



y Nx + yrj\ J 



+ 1^ + ^ Y^^ ]^h-'^ab 



+ 0{h^). 



P- =2x, 


(32a) 


Q- = -x, 


(32b) 


''' = ---2N' 


(32c) 


'•- = ---2N- 


(32d) 


A.-=. + |. 


(32e) 



Comparing this equation with the most general form of the linearized geometry tensor ([9]) and 
writing the parameter matrices in the form (jlOp , we read off the parameter values 

P+ = i2-N)x + y, 

, , N -2 y 

2 2iV' 

, 2-A^ 2-A^ 

These values of course satisfy the gauge- invariance conditions ()12p , since they result from a diffeo- 
morphism invariant gravitational action. 

Now we are in the position to follow the steps detailed in section IIIIBI to compute the PPN 
parameters of our theory. 

B. PPN parameters 

With the parameter values of the linearized theory obtained in ()32p we now compute the PPN 

parameters a^, 7^, 9^,a_^_ of our repulsive gravity model defined by ([26]) and ([27|) . The procedure 

for this is based on solving the equations of motion, for which we developed the technology in 

section IIIIBI A summary of the necessary steps is given at the end of that section. Choosing the 

standard PPN gauge such that 9~^ = yields the PPN parameters 

,1/3 4iV-4\ _ 1 7Nx + y 

a = ——- —- — , a 



3N \Nx + y Nx-y 
1 / 3 2iV-2 



?,N \Nx + y Nx-y^ 
+ _ (2JV - 4)x + 2y 



3NN^x^- 


y2' 










1 / 3 




iV- 
Nx- 


-2 

-y 


+ 


N 


3iV \Nx + 


y 


y 


1 


1 


5 








3Nx - 3y 




Ax 
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We now focus on the Newtonian limit of our theory. Recall that we demanded in [l|] a Newtonian 

limit where the gravitational interaction within each sector is attractive, while it is repulsive of 

equal strength between matter belonging to different sectors. This limit corresponds to the PPN 

parameters a"*" = 1 and a~ = — 1. These are achieved for parameter values 

2iV-l -2iV + 7 

y = wo .n • (34) 



6A^(2-iV)' ^ 6(2 -A^) 
Note that this recovers the same values that can also be obtained from a purely Newtonian cal- 
culation [l]. Substituting these values into ()33p simplifies the results for the PPN parameters 
to 

a^ = 1 , a~ = — 1 , 

+ - J_ - - 3 J_ 1 

^ ~iV' '^ ~2A^-7"^iV^2' 

Comparison with the observed values 7"*" = 1 and a']^ = — 2 displayed in ([7]) immediately shows that 
these are satisfied only in the case A^ = 1, i.e., when there is only one metric and a corresponding 
copy of the standard model, in which case our theory reduces to Einstein gravity. This is a 
dissatisfactory result since our aim was the construction of experimentally feasible gravity theories 
for N >l. 

This result shows that our model requires modification in order to match experimental bounds 
from solar system experiments. In the following we will make such improvements that adapt the 
theory to the observed values of the PPN parameters. 

C. Improved PPN consistent model 

Since the PPN parameters calculated for our theory do not match the observed values, it is 
natural to ask whether the theory can be modified so to reproduce the correct values. We will now 
show that this is indeed possible. 

We will modify the action and then proceed in complete analogy to the previous sections. 
First, we will compute the field equations from a variation of the modified action. Second, we will 
expand the metric around a flat solution and derive the linearized field equations. Third, we will 
read off the values of the parameter matrices and employ the linearized multimetric PPN formalism 
constructed in section IIIIBI We will only give a brief sketch of this calculation here. 
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We start from the gravitational action (j26p and add the following term which is consistent with 
our assumptions (i)~(iv) of section |II] that restrict the multimetric theories we consider in this 
article: 

N 

Sg = y^^ d'^x^og'"' {z~S\~S'\j + u~S'i¥j) . (36) 

^ 1=1 -^ 

This term contains two new constant parameters z, u that will be determined by PPN consistency 

below. The above modification is not the only possibility to achieve experimental consistency of 

our model, as we will discuss in the conclusion. 

Here, we continue by computing the equations of motion by variation. These take the form 

displayed in equations ([28]) . but the curvature tensor K^^^ in (f29|) attains a correction term ^^^ 

that is rather involved. We then compute the linearized curvature tensor. In addition to the result 

obtained in (j3ip this gives 



Next, we read off the modified parameter matrices P., Q , R, M.i l!L- Using the notation introduced 
in section IIIIAl we obtain the following modified results as compared to (|32p : 

P+ = (2 - N)x + y, P- = 2x, (38a) 

Q+ = ^L-lx-^, Q- = -x, (38b) 

, JV-2 y N-1 y z 

^ =^"-2iV + ^Ar^' ^ =-"-2iV-2iV' ^''^) 

M+ = ^^^x-^ + ^^^., M- = -x-^-^, (38d) 

2 2iV 2A^ ' 2N 2N' ^ ^ 

,, 2-iV 2-N iV-1 , ,^ y z-u , , 

Finally, we follow the steps detailed in section IIIIBI to compute the PPN parameters 
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a^,7^,^^, cr+- In comparison to ([33]) these take the modified values 



1 / 3 4iV-4 2{N-l){u-3z) 



3N \Nx + y Nx-y 

_ _ 1 / TNx + y 2{u- 3z) 
" ~ 3iV ViV2a,2_y2 + ^ 

_^ _ 1 / 3 2N-2 2(iV-l)(n-3z) 

^ ~ 3iV \Nx + y~ Nx-y^ H 

1/^3 iV-2 N {AN - 2)u + (6 - 9N)z + 3A^j/ \ 



?,N \Nx + y Nx-y y 

^_ 1 1 4n + 3(y-3z 



9+ = 0, 



3Nx - 3y 3y 3H 

(27V - 4)x + 2y _ Ax 



+ N^x^-y^ ' + iV2x2-y2- 

for H = 3(y2 + 2:^) _ 2Nx{u — 3z) + 2uy. One can now choose the parameters x, y, z, u so that not 
only the Newtonian limit agrees with our repulsive gravity requirement a^ = —a~ = 1, but also 
the observed PPN parameter values 7+ = 1 and a^ = —2 are obtained. These results are achieved 
for parameter values 

1 4-A^ 4-iV 12 -37V 

y = -o n^' ^ = —^ Jl^' ^ = -^ — 7^- (40) 



8 - 4iV ' " 8 - 4iV ' 8 - 4iV ' 8 - 4iV 

The remaining PPN parameters then are determined to be 7" = —1, 0~ = and a'^ = 2. 

Now that the PPN parameters of the improved model are consistent as far as we can determine 
from the linearized multimetric PPN formalism, one may ask whether the remarkable cosmological 
features of the original theory presented in [l|, such as the accelerating late time expansion and 
the big bounce, are still present in the improved version. 

It can easily be seen that this is the case by noting that the modification (j36p of the action 
is quadratic in the connection difference tensors S . Consequently the additional terms in the 
equations of motion caused by this modification contain at least one connection difference tensor. 
The remains of these are also seen in the linearized term (|37|) . For our simple cosmological model 
we assumed that a cosmological version of the Copernican principle holds in the sense that equal 
amounts of each type of matter are distributed homogeneously in our universe, and thus we could 
argue that all metric tensors of the cosmological solution should be equal at very early and very late 
times. It then follows that the connection differences vanish and the earlier obtained cosmological 
equations of motion are unchanged under the modifications we presented in this section. 



16 



V. CONCLUSION 

This article continued our discussion of multimetric gravity theories in which the gravitational 
field is described by N metrics, and which contain a corresponding number of standard model 
copies [l|, y]. The motivation to study these theories comes from the fact that they can be con- 
structed so that the astronomy of A^— 1 matter sectors appears dark for any given observer, and (for 
A^ > 2) so that different types of matter feature a mutual gravitational repulsion in the Newtonian 
limit. In this way they may naturally model cosmological effects such as late-time acceleration. 

From the particle theorist's point of view it can be regarded as a strength of our models that 
they do not introduce matter fields of unknown masses, charges or couplings. The non-gravitational 
particle content of the well-understood standard model is simply copied into the different sectors. 
Interactions between the sectors are mediated only through the coupling of the different metrics, 
as becomes clear from the action structure ([1]); the relevant cross sections will involve the New- 
ton's constant squared. Hence, direct experimental observation of the other matter types will be 
extremely difficult. 

To discuss the gravitational field content of our models, we repeat an observation from [4]: while 
N symmetric two-tensors /i^^ appear on the linearized level, diffeomorphism-invariance merely 
implies a single gauge symmetry under S^h^^f^ = 2d(^a^i,^ ^or common gauge parameters^,- This 



231]. Since 



type of gauge symmetry is required for the definition of a massless particle of spin two 

every observer in our theories can choose to relate the gauge symmetry to his own metric field, 

we may interpret our theories as containing one graviton and further A^ — 1 symmetric two tensor 



m 



fields that cannot be interpreted as spin two particles p 

In this article we asked whether theories of this type are consistent with data available from 
high precision experiments at the solar system level. A suitable framework to test gravity theories 
based on a single metric exists, and is known as the parametrized post-Newtonian formalism. In 
this article we have constructed a simple extension of this PPN formalism to multimetric gravity 
theories. We worked on the assumption that regions exist which are dominated by matter belonging 
to a single standard model copy, and where one may neglect the influence of matter from different 
sectors on the gravitational field. This is a reasonable assumption for repulsive gravity theories of 
our type because different types of matter should separate as the universe evolves. 

Our multimetric extension of the PPN formalism features an additional set of PPN parameters 
describing the mutual gravitational interaction of matter belonging to different sectors. We have 
shown that a subset of these parameters can be obtained already from the linearized equations 
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of motion. This yields a quick method of selecting consistent theories or discarding unphysical 
ones. We have applied the linearized multimetric PPN formalism to the particular gravity theory 
presented in p|. We computed the PPN parameters and found that these do not agree with 
the observed values. We finally could present an improved theory which now does agree with 
observations while still featuring the same promising cosmological effects as the original theory, 
most notably the accelerating late-time expansion of the universe. 

It is worth pointing out that the particular improvement (j36p of our multimetric repulsive gravity 
model is not unique. Two new parameters control the modification terms and need to be determined 
by experimental consistency requirements. Even within the class of quadratic connection difference 
terms, one could discuss other examples of modification, with even more additional parameters, 
that also yield correct PPN values without changing the cosmological properties. The linearized 
multimetric PPN formalism developed in this article does not well distinguish theories of this type. 

Hence it would be desirable to find and establish further physical and mathematical principles 
in order to restrict the possible terms in multimetric gravity actions. The most obvious restriction 
one should pose is that the remaining PPN parameters that cannot be obtained from the linearized 
theory should also agree with the observed values. The determination of these parameters would 
require a perturbative expansion of the equations of motion up to quadratic order in the metric 
perturbations which is very involved for multimetric gravity theories, but should be carried out in 
future work. One mathematical idea to restrict possible gravitational actions could be to enlarge the 
symmetry group. While we restricted to a discrete exchange symmetry with respect to arbitrary 
permutations of the sectors, one could think of establishing a continuous symmetry group that 
mixes the sectors. In any case, the search for further principles behind multimetric gravity actions 
is still an open question. 

Appendix A: CoefRcients of the linearized PPN ansatz 

In this appendix we display the detailed expression for the coefficients c^ , . . . , Cg used in the ex- 
pansion of the geometry tensor given in equation ()18p . These can be computed using the expression 
for the linearized geometry tensor ([9]) and the linearized PPN metric ansatz ^ 



4 = -{N+ + Q+)a+ -{N- 1){N- + Q-)a- 

+ (M+ + 3iV+)7+ + {N - 1){M- + 3iV")7- (Al) 

- (M+ - N+)e+ -{N - i){M- - N-)e- , 
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q = -(A^^ + Q-)a+ - (iV+ + Q+ + iN- 2){N- + Q-))a^ 

+ {M- + 3N-)-f+ + (M+ + 3iV+ + {N- 2){M- + 3N-))-f- 
- {M- - N-)6+ - (M+ - N+ + {N- 2){M- - N-))e- , 

c+ = Q+(7+ + (iV-l)Q-c7;, 
c- = Q'al + (Q+ + (iV - 2)Q-)a^ , 

2c+ = (P+ + 2i?+)a+ + (A^ - 1)(P" + 2R-)a- 
- {P+ + 6i?+)7+ - (iV - 1)(P- + 6R-)-f- 
+ (P+ - 2i?+)0+ + {N - 1){P- - 2R-)9- 
+ (P+ + 2Q+)cj+ + (Af - 1)(P' + 2Q~)aZ , 



(A2) 

(A3) 
(A4) 



(A5) 



2c^ = (P- + 2R-)a-+ + {P+ + 2R+ + {N - 2){p- + 2R-))a- 
- {P- + 6i?")7+ - (P+ + 6P+ + (iV - 2)(P- + 6P"))7" 
+ (P~ - 2R-)e+ + (P+ - 2P+ + {N - 2){P- - 2R-))0- 
+ (P- + 2g-)cj+ + (P+ + 2Q+ + iN- 2){P- + 2Q-))a-_ 



(A6) 



= iV+a+ + (iV- l)iV"a~ 

- (M+ + 3Af+ + Q+)7+ -{N- 1)(M- + SiV" + Q-)r 
+ (M+ - iV+ - Q+)^+ + {N - 1){M- -N- - Q-)e- , 



(A7) 



N-a+ + (iV+ + (iV - 2)N-) 



a 



- (M- + 3N- + Q-)7+ - (M+ + 3iV+ + Q+ + (iV - 2)(M- + 3N~ + g~))7' 
+ (M" -N- - Q~)e+ + (M+ - Af+ - g+ + (Af - 2){M- - N' - Q"))0" , 

c+ = P+a+ + (A^ - l)p-a- 



- (P+ + 3P+)7+ - (Af - 1){P- + 3P~)7" 
+ (P+ + 2Q+ - P+)^+ + (AT - 1)(P- + 2Q- 



P- 



(A8) 



(A9) 



= P~a+ + (P+ + {N - 2)R-)a- 

- {P- + 3P")7+ - (P+ + 3P+ + {N - 2){P- + 3P"))7" 
+ (P- + 2Q- - p-)^+ + (P+ + 2Q+ - P+ + (Af - 2)(P- + 2Q- - p-))0- . 



(AlO) 
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